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We compute the inclusive cross-section of fa tensor mesons production in proton- 
CN ' proton collisions at high-energy. We use an effective theory inspired from the tensor 

meson dominance hypothesis that couples gluons to fa mesons. We compute the 
differential cross-section in the fcj_-factorization and in the Color Glass Condensate 
^| formalism in the low density regime. We show that the two formalisms are equivalent 

i—i, for this specific observable. Finally, we study the phenomenology of fa mesons by 

Q H i comparing theoretical predictions of different parameterizations of the unintegrated 

Q-T gluon distribution function. We find that /2-meson production is another observable 

& '. 

. that can be used to put constraints on these distributions. 
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0\ . I. INTRODUCTION 

-xf ; 

High energy hadronic collisions are complex phenomena that combine many-particles 
\ physics and Quantum Chromodynamics (QCD). Making predictions for the production of 
particles in these reactions involves the understanding of both hadron wave-function and pro- 

> ■ 

cesses of particle creation. Many advances have been made in the last decades in those fields 
through the application of perturbative QCD (pQCD) to the description of experimental 
data. The pQCD analysis relies generally on factorization schemes like Collinear Factor- 
ization and /c^-factorization. Different procedures are implemented in these approaches to 
improve the perturbation expansion by resumming infrared divergences. Many observables 
are well-described by QCD factorizations. For example, /c^-factorization was used success- 
fully in proton-antiproton collisions for heavy-quarks production 

while collinear factorization is one of the main computational tool for deep inelastic scat- 



tering and for a number of other applications 10l lllj. The main difference between the two 



formalisms is their validity range, /c^-factorization describes semihard processes character- 
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ized by typical momentum transfer \x obeying A.q CD C /i 2 C s while hard processes with 
A.q CD <C fi 2 ~ s is the domain of collinear factorization. 

Some time ago, new ideas were developed to take into account effects due to the recom- 
bination of gluons in nuclei at very high energy (or small-x) 12, 0, Q, 0, Q- These 
saturation effects introduce a new scale called the saturation scale Q s at which the proba- 
bility of gluon recombination becomes important. A naive estimation of Q s shows that it 
depends on the momentum fraction x and the number of nucleon like Q 2 ~ A 5 x~ x 0, fl5 ]. 
At small enough x or large enough A, the saturation scale becomes hard {Q 2 S ^> Aq CD ) 
and weak coupling techniques can be used. The Color Glass Condensate (CGC) includes 
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3,0, HQ- 

More precisely, it 



these saturation effects in a semi-classical formalism 
resums the large order contributions coming from the interaction of gluon cascades in the 
regime where A.q CD <C /x 2 < Q 2 <C s. In proton-proton (pp) collisions at RHIC energy (with 
Q s < /i), saturation can be neglected and collinear or fc^-factorization can be used in the 
calculation of observables. It is interesting to compare the predictions of CGC and factor- 
izations in that regime. For heavy quark production [l7\ and gluon production 18, 19], the 



CGC in the low density limit and fc^-factorization are equivalent in pp collisions. We show 
in this paper that this also hold true for tensor meson production. 

In this study, we focus on the production of /2(1270), the lightest spin-2 meson. It is a 
bound state of light quarks (u or d) and antiquarks, it has a mass of 1.27 GeV, a decay width 



of 185 MeV and quantum numbers of I G (J PC ) = + (2 ++ ) 20]. This partic 



e was studied in 



the past using different methods such as chiral perturbation theory [211 . [22], [23j and effective 
field theory [24| . An other approach is to use the tensor meson dominance (TMD) hypothesis 
where one assumes that the energy-momentum tensor form factor are dominated by the 
exchange of an f 2 meson. The idea of TMD was first applied to phenomenology by Renner 



25l . |26| to describe the tensor-meson exchange channel in pion-nucleon scattering. This idea 
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33, 



34, 



35| . Borrowing ideas 



was then used in many other applications 
from the TMD hypothesis, namely the coupling of spin-2 mesons to the energy-momentum 
tensor of strongly interacting matter (including quarks, gluons and hadrons), we show how 
the inclusive cross-section of f 2 mesons can be computed in proton- pro ton collision at high- 



energy (RHIC and LHC) . This type of coupling is used in 34| and in 36] for the computation 
of the f 2 decay rate into photons and pions. In the study described in 36(, the authors 
compute successfully (within experimental error) the decay rate of f 2 into two photons 
using ADS/QCD. To obtain this accurate result, they must include both the gluon and quark 
contributions to the energy-momentum tensor. In this paper, we consider only the gluonic 
part of the energy-momentum tensor since the hadrons in high-energy collisions interact 
mostly through gluons. More precisely, we study the process where the interaction of two 
gluons coming from different nuclei interact and produce one on-shell f% meson because this 
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is the dominant process at high-energy. 

We consider only the case of proton-proton (pp) collisions and discuss briefly the case 
of proton-nuclei (pA) collisions. For nuclei-nuclei (AA) collisions, the total number of Ja- 
mesons produced should be important but most of them cannot be detected because they 
decay inside the medium created by the collision. This can be seen as follow. Our analysis 
shows that tensor mesons are created during the first instants of the collision (t < 1 fm/c), 
when the system is still in non-equilibrium and well described by strong classical field. RHIC 
data indicates that from 1 fm/c up to 10 fm/c, a thermalized medium is created. The mean 
lifetime of f 2 mesons is about 1.1 fm/c, which is smaller than the lifetime of the medium. 
Moreover, the probability that these particles travel distances larger than the radius of the 
nuclei (or larger than the medium) is very small unless their momenta obey |p| 3> M. This 
is very unlikely because the typical transverse momenta will be M ~ 1.27 GeV and the f 2 
are created by small-x gluons which have very small longitudinal momenta to begin with. 
Thus, the f 2 will mostly decay inside the medium and the decay products will rescatter, 
losing information about their origin. In that sense, most tensor mesons produced should 
be unobservables in AA-collisions unless they are produced on the surface. This can be 
seen experimentally in the invariant mass distribution where the f 2 signal is essentially non- 
existent 37] . This is not the case for pp and pA because no medium is created and the 



particles obtained from f 2 decay can escape without rescattering. For these reasons, our 
present analysis can only be applied to pp and pA collisions. 

This paper is organized as follow. Section 2 is dedicated to the statement of the spin- 
2 tensor mesons effective theory. We use an effective free Lagrangian describing spin-2 
particles dynamics that is similar to Kaluza-Klein modes in extra- dimensions studies. The 
f 2 interacts with strongly interacting matter through the energy-momentum tensor. In 
section 3, the computation of the cross-section at leading order in perturbation theory using 
fc_i_-factorization is performed. The process considered is the tree-level interaction of two 
gluons giving one f 2 . In section 4, the same calculation is done using the CGC. We start by 
deriving a reduction formula relating the cross-section production to a correlator of energy- 
momentum tensors that can be evaluated using the CGC formalism. We use the solution of 
the gauge field in covariant gauge at leading order in the color charge densities to evaluate 
the cross-section. The two formalisms are then compared and shown to be equivalent in the 
limit of low densities (dilute systems). In both cases, the production cross-section is related 
to the unintegrated gluon distribution functions of the two protons. Section 5 is devoted 
to the phenomenology of f 2 meson production at RHIC. In that section, the predictions of 
standard parameterizations of the unintegrated distribution function are compared against 
each other. 
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Throughout the paper, we use both light-cone coordinates defined by 

, p° + p 3 _ p° — p 3 

p = — ^— ; p = — ^— (i) 

y/2 V2 

and Minkowski coordinates. It should be clear by the context which one is used. We also 
use the metric convention rj^ = (1,-1,-1,-1). 



II. EFFECTIVE THEORY 

In our approach, the tensor mesons couples to the rest of the strongly interacting matter 
through the energy-momentum tensor. The energy-momentum tensor is a conserved current 
that has the same quantum numbers as tensor mesons, so a coupling of this type is natural. 
The effective interaction Lagrangian descri bing the interaction of gluons, quarks and hadrons 



with tensor mesons is simply given by 34. 136] 



Ant(^) — fnv( X )^uovs( X ) + ,,. /^( X )^quarks( X ) + fl»> ( X ) ^hadrons ( X ) (^) 

where f^ u is the symmetric spin-2 tensor field, k and are the coupling constants and 
^giuons,quarks,hadrons are the energy- momentum tensors of gluons, quarks and hadrons respec- 
tively. This type of interaction preserves all the symmetries of the matter Lagrangian by 
construction. For simplicity, we consider only one species of tensor meson, namely /2(1270). 
Generalization to include other species like the /2(1525) is straightforward. Note also that in 
this theory, the coupling constants k are free dimensionful (dimension of energy) parameters 
that need to be fixed by experiments. It is possible to fix their value by assuming that they 
are the same (k = K q = k^). Within this assumption, the hadronic and partonic sectors of 
the theory couple to tensor meson in the same way. Then, k can be fixed using the f% pions 



decay and one finds that k ~ 0.1 GeV 30j, |34j, [36|. The theory used in this paper differs 



from this approach, we use different coupling constants for each sectors Kh). The 

value of the coupling is then evaluated by comparing to experiments (see the last section for 
a comparison with STAR data). 

For a physically consistent formulation of f% tensor Lagrangian, we borrow from recent 
developments of Kaluza-Klein theory and write the following free Lagrangian describing 



spin-2 particles dynamics [38l. l39l| . 



Uee(x) = -lf^( X )(d? + M 2 )f„ v (x)+~r p (x)(d 2 +M 2 )fZ(x) 

- r(x)d,d l/ r p (x) + r(x)d fl d a f^ X ) (3) 

The first term of the Lagrangian is the usual kinetic term while the other terms are necessary 
to have the right number of degrees of freedom. This will be discussed in more details in 
section IIV Al 
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This non-renormalizable Lagrangian describes a low energy effective theory for tensor 
mesons dynamics. Therefore, its domain of validity is restricted to a certain energy domain 
which is of the order of the mass of f'2 mesons because it is the only scale appearing in the 
problem. Although this cannot be proven rigorously, the results obtained from this theory 
should be taken with caution at very high transverse momenta. 

In high-energy proton-proton collisions, it is well-known from parton distribution func- 
tions that gluons dominate the cross-section, so the main contribution to f'2 production 
comes from the gluonic sector of the theory. The energy-momentum tensor of gluons can be 
obtained from the Yang-Mills Lagrangian by varying the metric and is given by 

TSL.{z) = ^G apa (x)G:"(x) - G£{x)G?(x) (4) 

where 

G?(x) = - PA&x) + gf abc A%(x)A»(x) (5) 

is the usual field-strength tensor and is the gauge field of gluons with color index a. The 
Feynman rules associated with this Lagrangian are presented in Appendix [A] This effective 
theory can now be used to compute the cross-section of production in high-energy pp 
collisions. 



III. /2-MESON PRODUCTION CROSS-SECTION FROM PQCD 

In this section, the / 2 - m eson cross-section in proton-proton collisions is evaluated in the 
fc^-factorization formalism. This formalism can be used when the collision is semihard, 
meaning that Aq CD <C /i 2 <C s where s is the squared center of mass energy, Aqcd ~ 
200 MeV is the usual QCD scale and \i is the typical parton interaction scale. The parton 
interaction scale is related to the transverse mass of the produced particles fi 2 ~ M]_ = M 2 + 
k\. For relatively low transverse momenta as the ones measured in experiments (for example, 
for momenta measured up to < 10 GeV), the interaction scale for production is 
within M 2 < ji 2 < 100 GeV 2 , which obey the inequality for semihard processes. Clearly, 
the fi observed at RHIC and LHC should be produced from these semihard processes, 
unless they are measured at transverse momenta of the order of the center of mass energy 
(\k±\ ~ a/s). Moreover, these processes occur at small-x (x ~ M]_/s 1) where the collinear 
factorized perturbation theory is spoiled because of logarithmic terms like [ln(/x 2 /Ag CD )a s ] n , 
\\n([i 2 / Aq CD ) ln(l/o;)a s ] n and [ln(l/x)a s ] n 0, Q, 0, Q] . The essence of fcj_-factorization is 
to resum these "large-log", leading to a description in terms of the unintegrated gluon 
distribution functions <f>(x, k±, /z 2 ). These functions give the probability tofind a gluon with 
longitudinal momentum x and transverse momentum k± at the scale fi 2 



HQ. 
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It should be noted that fc^-factorization have not been proven rigorously for $2 production. 
However, it is analogous to proven heavy quark production described in and to Hide's 

production studied in [g]. Moreover, as shown in section HVl the Color Glass Condensate can 
be used to justify this approach since the two formalisms give the same cross-section. This is 



similar to the work described in [17| where it is shown how to recover fc^-factorization from 
the CGC formalism. For these reasons, using /c^-factorization for /2-production at RHIC 
and LHC energies is justified. 

A. Cross-Section in fc^-factorization 

The inclusive cross-section in the /c_|_-factorization formalism is given by 0, Q, S] 

(n ,3 OT? da pp ^ f2X 2 f 1 dx 1 dx 2 f d 2 q±d 2 p± 2 2 2 2 

[2^Y2E k — — = 16vr / / 0i(xi,p ± ,/x )<f> 2 {x 2 ,q ± ,ii ) 

a k Jq x\ X2 J 

*(2*) 3 2E k d3k (6) 

where 

(2vr) 3 2E fc — = -\M^^\ 2 W5\p + q-k) (7) 

is the high-energy limit of the cross-section for off-shell gluons g* to on-shell f2, s = X1X2S 
is the fc^-factorization flux factor X\2 are momentum fractions of gluons and 

01,2(^1,2, k±, fi 2 ) are unintegrated gluon distribution functions of proton 1 and 2. The unin- 
tegrated distribution functions are related to usual parton distribution functions of gluons 
(appearing in collinear factorization ) by 

/ dk 2 ± <p(x,kl,n 2 ) « xG(x,fi 2 ) (8) 
Jo 

where G(x,fi 2 ) is the usual gluon distribution function in collinear factorization. 

B. Kinematics and Matrix Element 

To compute the production cross-section of /2- m esons, the high energy limit of the lowest 
order matrix element between two off-shell gluons and one on-shell f 2 has to be calculated. 
The Feynman diagram included is shown in Fig. (DO). It can be evaluated by using the 
Feynman rules presented in Appendix |A] The main difference with usual pQCD calculation 
is in the average on polarization of off-shell gluons. In /^-factorization, the polarization 
tensor obeys [J 0] 

J2^(pyx(p) = ^ (9) 
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Parton 1 



f 2 

k 2 



Parton 2 



FIG. 1: Feynman diagram included in the lowest order calculation of f 2 production in the k±- 
factorization formalism. Two partons with momenta P and Q, each part of proton 1 and 2, emit 
two off-shell gluons with momenta p and q. The gluons combine to produce a ji meson with 
momentum k. 

where p± = (0,p±, 0). The sum on polarizations differs from the usual result because we are 
considering off-shell gluons with a virtuality given by p 2 = —p\. The exact form is due to 
the coupling of gluons to partons through eikonal vertices as well as gauge invariance and 
Ward identities [2]. This can also be shown to be gauge invariant 

In the center of mass frame, the 4-momenta of partons inside the proton moving in the 
±z direction in Minkowski coordinates can be written, as: 

*=(£o,o,f);g=(£co,-£) (10, 

Then, the momenta of gluons in the large energy limit (\p±\, \q±\ *C is simply 

v = {—>p±>—) ; *={—>k>-—) ( n ) 

The / 2 - m eson is on-shell and obeys k 2 = M 2 . 
The matrix element for f 2 production is 



IU9V-/2I2 = 1 V V" |7-ffV^/2 12 



(N 2 

K c 'ah A,A',A" 



„, (iV?-l)2— pi ql 

xV^(k,q,p)[v^(k,q,p)]* (12) 

where a, b are color indices, N c is the number of color, V£f vl (k,q,p) are vertices defined in 
Appendix [A] and the projection operator is defined as 



A 

with g^ = -k^k v /M 2 . 



(13) 

k 
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By contracting the tensors p± p ,p± v and g± C r,5 , ± 7 with the vertices, the expression of the 
matrix element can be simplified further (with p + and q~ in light-cone coordinates): 

w~*r - w^^ ^'^J^ (14) 

where we defined H t f u (p±, q±) as 

H]_\p ± , q ± ) = -Hl 2 (p ± , q ± ) = q'p 1 - q 2 p 2 (15) 

H[ 2 (p ± , q ± ) = Hf(p ± , q ± ) = q'p 2 + q 2 p l (16) 

All the other components of are zero. 



C. Cross-Section 

Once the matrix element is known, it is possible to evaluate the differential cross-section. 
Using the kinematics (Eq. (fTTj) ) and the delta functions, the integrations on x li2 can be 
easily done. The /c^-factorized differential cross-section is 

(o \s OJP d ° W ^ hX 1R 4 P^ap(k) f d 2 q ± d 2 p ± 2 2 

(2?r) 2Ek ^f— = 16?r (N? - IW J ^]^ (27 ° 6 iP± + q± ~ k±) 

xfr ^, ^)0 2 (x_, gi, ^ HTiP± > q±) 2 H 2 fiP± > q±) (17) 

where x± = [E^ ± fcJ . This is one of the more important result of this paper. It relates 
the production cross-section of / 2 mesons to unintegrated gluon distribution functions. The 
phenomenology of this equation and its derivation from the CGC formalism are studied in 
the next sections. 



D. Limit of Collinear Factorization 

The procedure to recover collinear factorization cross-sections from /c^-factorization is 
well-known jl, 0, 0, [l7] and will serve as a consistency check for Eq. 0171) . The limit 
\p±\, |<Ll | — > has to be taken in the matrix elements and the integration on the azimuthal 
angle has to be performed. The last step is to use the relation Eq. (jSj) to make the last 
integral and relate the unintegrated distributions to the collinear distributions. 

For the case of / 2 production, the integration on azimuthal angles of the p±, g^-dependent 
part of the matrix element is 

27T d9 p d9 k lim P^(k) H ^±^)n?(P±><!±) =Sn * (18) 
o b±U<?±Ho P±Q± 
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Combining this result with Eq. (J5J) and Eq. (|T7|) . the cross-section becomes (s is the center 
of mass energy) 

{2n)32Ek ^¥r- = (nT-^s W^M**' f)G{*-^ 2 ) (19) 

which is exactly the same as the expression computed in the collinear formalism in Appendix 
|B] (see Eq. IB6j) . Thus, Eq. (1171) has the right collinear limit. 



IV. PRODUCTION OF / 2 -MESONS FROM THE CGC 

In collisions at very high energy, the wave function of nuclei is dominated by soft gluons 
( small- x where x is the momentum fraction). The Color Glass Condensate (CGC) is a 
formalism that describes the dynamics of these degrees of freedom in the limit of large 
center of mass energy. In this approach, hard partons that carry most of the longitudinal 
momentum, and soft gluons that have small longitudinal momentum component, are treated 
differently. Hard partons act as sources for soft gluons and are no longer dynamical degrees of 
freedom. The occupation number of soft gluons is large because of an emission enhancement 
at small-x, so classical field equations can be used to understand their dynamics (for reviews 
of CGC, see 0, Q, Q ) ■ 

In the CGC, computing a physical quantity involves two main steps. The first one is to 
solve the Yang-Mills equation of motion 

[D fl ,F^(x)] = r(x) (20) 

where the current J% (x) = d u+ 8(x~)pi n (x \ ) + 8 u ~5(x + )p 2 ^ a (x±) represents random static 



sources localized on the light-cone [14|, 115]. In that case, P i, 2 (x±) are color charge densities 
in the transverse plane of proton 1 and 2 respectively. The next step is to take the average 
over the distribution of color charge densities in the nuclei with weight functionals Wi^Pi^] 
that include the dynamics of the sources. For any operator that can be related to color 
charge densities, this can be written as 

(0) = Jvp 1 Vp 2 0[p u p 2 }W 1 [x uPl }W 2 [x 2 ,p 2 } (21) 

Computing the weight functional is a highly non-perturbative procedure so it usually involves 
approximation based on physical modelling. In the limit of a large nuclei at not too small x, 
it can be approximated by the McLerran-VenuKo pal an ( MV) model, which assumes that the 
partons are independent sources of color charge [12j, [l3|. In pp collision, this kind of model 
fails and a more phenomenological approach, through unintegrated distribution function, is 
better suited. 
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One important ingredient is missing for the computation of /2-meson production cross- 
section. It is necessary to have a relation between the cross-section and a correlator that 
can be evaluated using Eq. (j2T]l . This is done in the next section by deriving a reduction 
formula for our effective theory. 



A. Reduction Formula for the Effective Theory 

The production of tensor mesons from the CGC can be calculated from a reduction 
formula derived from the Lagrangian Eq. ([3]). The first step is to compute the equation 
of motion of the tensor field f^ v given as usually by L^ u (x) = Sf ^ (x) = 0. Using L^ u as 



well as L^, d^L^ and d il d v L VLV (see [38|] for details), the following equation of motion is 
obtained (note that for the simplicity of this derivation, we chose only the coupling to one 
type of particles so that the interaction Lagrangian looks like \f^ u (x)T^(x), the case with 
the Lagrangian given in Eq. ([2]) can be easily dealt with at the end): 

(d 2 + M 2 )f^(x) = -e^(x) (22) 



where 



O^x) = T, v {x) - + g, u ) ^P- (23) 



as well as the following two constraints: 



d »f {x ) = -MM and Mx) = (24) 
H y ' 3M 2 k Jpy ' 3M 2 k v ; 

These constraints make f^ u to have the right free-field limit with only 5 propagating modes. 
The source term satisfies 



p 

and 



©^) = -^(5 2 + m 2 )t;(x) (25) 



d^(x) = -^(d 2 + M 2 )d u T^x) (26) 

The above equations can also be interpreted as operator equations. Since the equation of 
motion is linear in f^, the general solution of the equation of motion Eq. ([22]) is easy to 
obtain. The Fourier transform of the space coordinates to momentum-space is taken and 
the solution for the field operator is given by 

/^(t,k) = /W(t,k) + - / dt'G rct (t-t',k)Q, u (t',k) (27) 
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where fjh)(t, x) satisfies the free field equation and the retarded Green function in the mixed 
t and k representation is 



GU(t , k ) = !E^fl 9 ( t ) 



(28) 



It can be easily shown in momentum space that the solution of obeys the constraints in 
Eq. (1241 . Therefore, this solution obeys the requirements of spin-2 particles and describes 
the right dynamics. 

The inclusive cross-section is related to the average number n of ji produced. In terms 
of creation and annihilation operators, n is given by 



dn 



(27r) 3 22^ = £<iiiit| 



*"]<; JotlH U 'k/OUt 



|init) 



(29) 



where A = 0, ±1,±2 are the polarization of the spin-2 tensor field. Here |init) specifies 
a quantum state at t = — oo and the subscript 'out' specifies that the operator is to be 
evaluated at t = oo. As t — > oo, we expect that the density of particles will become low 
and the fields become asymptotic. We start by writing the creation/annhilation operators 
in term of the field like 



t— >oo 



a£) out = lim ((£,(*))* E k r{t,k)+id t r(t,k) 



where the polarization tensor is traceless, transverse and satisfies 
Using the solution Eq. (|27j) . we then get 



(30) 



(31) 



(&k) 



,, ut = lim(a£e iEkt ) ilcc 

t— >oo 

+ hm- f dt'e-^^{e^(k)YQ, u (t',k) 



where we used Eq. ( 128|) to get 

iG Iet (t,k) + E k G rct (t,k) = ie- iEkt 9(t) 

Here. 



{aie~ iEkt ) 



free 



(k))* \E k f^ v (t,k)+idJ^(t,k) 



(32) 



(33) 



(34) 



describes the time evolution of the initial tensor meson population. Eq. (I3"2l of course is a 
reduction formula for f 2 production. To simplify the analysis, it is assumed that the initial 
state contains no or little tensor mesons 



(a^frce |init) = 



(35) 
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which is not an unreasonable assumption for high energy hadronic collisions. The average 
number of fi produced is then simply 

{27r)32Ek d^k = pW (*)< e U*)<W*)> (36) 

where Q pa (k) is the Fourier transform of Q pa (x) evaluated at a tensor meson on-shell mo- 
mentum. The projection operator P^ vpa {k) is obtained by summing over polarizations and 
is given by Eq. (|T3|) . The angular brackets (O) here indicates expectation value of O in the 
initial state. Due to the presence of the projection operator in Eq. (I3T)|) . some more simplifi- 
cations can be made. The partial derivatives in Q pu simply results in factors of k p k v in the 
Fourier transform. This implies that in the expression (1361) . those partial derivative terms do 
not contribute due to the transversality of the projection operator {P^ upa {k)k ll = 0) when 
k is an on-shell momentum. Furthermore, since P^p" i s made up of eig, it is also traceless 
in the sense that P^ v p° 'g^ v = P^ 9 " g pa = 0. Hence the g^ v term in does not contribute 
to this momentum distribution, either. By making these simplifications and by doing the 
same derivation with the interaction Lagrangian in Eq. (T5]), the following simple result is 
obtained: 

W2E k ^- = ^P^mfl v {k)T hP M) (37) 

i,j J 

where the indices are i, j = {gluons, quarks, hadrons}. This is the main result of this section. 
It relates the average number of mesons produced to a correlator of energy-momentum 
tensors. This correlator can then be evaluated using any analytical or numerical meth- 
ods available. The average (...) depends on the system studied. Looking at a plasma at 
equilibrium, it could be computed using finite temperature field theory or even ADS/CFT. 
However, as stated in the introduction, tensor mesons cannot be observed if the thermalized 
medium lifetime is longer than the decay width of f2 mesons. Therefore, these techniques 
are not pursued here. In high energy hadronic collision, the correlator can be evaluated 
using the CGC once T pv is expressed in terms of the sources p\^- This is done in the next 
section. 

The result obtained so far is exact given the Lagrangian ([3]). The only assumptions made 
along the way are: 

1. There is no tensor mesons in the initial state 

2. The tensor meson content of T^ u is ignored (T^ ~ 0) 

3. The tensor meson can be produced on mass shell 
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The first assumption is justified by the fact that in high-energy collisions, the number of 
tensor mesons in a nuclei before the collision is negligible. The second assumption is more 
subtle as it does not appear explicitly in this calculation. If this assumption is not made, 
the energy-momentum tensor depends on f^ u through the hadronic sector so the equation of 
motion becomes non-linear and much harder to solve. In a high-energy collision, the cross- 
section is dominated by gluons, so it is consistent to assume that -C T^ ons as long as the 
/2 mesons are produced in a small number. Finally, the last assumptions is that tensor meson 
can be produced on-shell so that their spectral density is p(M 2 ) ~ 5(p 2 — M 2 ). However, fa- 
mesons are resonances, so the spectral density should look rather as a Breit-Wigner function 
p(M 2 ) ~ T/((p 2 - M 2 ) 2 + M 2 r 2 ) where T is the decay width. Since M 2 > T 2 , we expect 
that corrections due to the finite decay width can be neglected and that we can approximate 
p(M 2 ) ~ T/((p 2 - M 2 ) 2 + M 2 r 2 ) « 5(p 2 - M 2 ). 

Having expressed the average number of /a produced in terms of a correlator of energy- 
momentum tensors, it is possible to compute the cross-section in the CGC formalism. In 
high-energy collisions, we have that T^ ons > ^quarks' ^Tadrons so we can ne gl ect the quark 
and hadron contributions. Thus, in the CGC, the inclusive cross-section for f 2 production 
is given by [l^, 0] 

da _ f J2u fn-\3n^ dfi ( b ±) 



(2ir) 3 2E k ^- = j d 2 b ± (27r) 3 2E, 



d 3 k 

-^P^upa(k) J d 2 b ± J V Pl V P2 f ^ ns [Pi > P^gLns [ Pl , p 2 ] 

xWMWtfaM (38) 



In this expression, b± is the impact parameter. The energy-momentum tensor TgJ uons is a 
functional of the source once the Yang-Mills equation of motion of the gauge field is solved 
(see Eq. (j4]) and (jSJ) for the expression of the energy-momentum tensor as a function of the 
gauge field). 

B. Calculation of the Cross-Section and relation to fc^-factorization 

In this section, Eq. (!38j) is evaluated explicitly to leading order for pp collisions. The 
expansion parameters in a dilute system like a proton are the weak color charge densities 
obeying pi^/^l *C 1. The first step is to solve the Yang- Mills equation. In a full solution, the 
gauge field would be a functional to all orders in the color charge densities. However, because 



the sources are weak (or equivalently, if the typical transverse momentum is large [ly, |46| ) . 
the solution can be truncated and the calculation can be done analytically. The solution of 
Yang-Mills equation to first order in pi^/k 2 ^ was found in covariant gauge (d^A^x) = 0) in 
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HQ. It is given by A»(k) = A$(k) + A£(k) + A^k) where A${k), A%(k) are 0( Pl ) and 
0(p 2 ) respectively, but A^ 2 (k) is 0(pip 2 ). The explicit components of the gauge field are 

k 



AUk) = 2n6(k~)^> (39) 



A 2>a (k) = 2n5(k^)^^ (40) 

All the other components are zero. The field A± 2 gives higher order contribution to the 
energy-momentum tensor correlator (T^ v \k)T pa (k)) , so it is not considered here. The lead- 
ing order contribution to the energy-momentum tensor correlator is 0(pfp 2 ) and that can 
only come from different combinations of A± and A 2 . A term containing A± 2 is at least 
O(plpl) or 0{p\p\). In principle, terms of O(pf), 0(p 2 ), 0(pip 2 ) and 0(pfp 2 ) can also be 
included but it can be shown that they do not contribute to the cross-section. Finally, the 
non-Abelian part of the energy-momentum tensor can also be neglected because it gives 
contributions of at least 0(p\p%) or 0(p\p\). Basically, the leading contribution of the 
energy-momentum tensor correlator in a dilute system is given by the Abelian part of the 
energy-momentum tensor computed with the gauge field A^ eading = A± + A 2 . 
The Abelian part of the energy-momentum tensor in momentum space is 

+ Pa (k - p) a A»(p)A» a (k - P )- Pa (k- P yAZ(p)A a a (k - p) 

-p»(k - p)X(pK(* - P) + P"(* ~ pYA a , a (p)A a a (k -p)\ (41) 



According to our previous discussion, the correlator (T^ u (k)T pa (k)) at leading order is given 
by changing A p — > Ay eading = A± + A 2 in the expression of the energy-momentum tensor. 
Then, the computation involves the evaluation of correlators with four gauge fields that have 
the following general form 

d j0f^(P, k, q)(AiUp)<a(k - P )Ai tb (q)A^ b (k - q)) (42) 

where = 1 or 2 and /^y^(p, k, q) is a function whose value is determined by the Lorentz 
indices structure and the momentum factors appearing in front of the fields inside the ex- 
pression of T p,u {k). We need to evaluate the correlator 

XZ: ^.,Ak.p. q) ee {A\l a {p)A\la{k ~ p)A p S3jb (q)A^ b (k - q)) (43) 

for all values of s< using the expressions of the gauge field in Eqs. (I3"9"j) and (T4"0"j) . They are 
given by 

A1%i,2(k,p,q) = 5 y - 5 P+ 5 a - {2ti f 5(p-)5(k + -p + )5(q~)5{k + -q + ) 

x (p{ a (p±)pi,b( ( l±)) (Pl a (k±-p±)p2,b{k± - q±)) 
p\q\ (k± - p±) 2 (k± - q±) 2 



14 



5> 1+ 5 u -5 fJ -5 r7+ (2ir) A 5(p-)5(k + - p + )5(q + )5(k~ - q 
(p\ ja (p±)pi,b(k± - q±)) (p{ a {k± -p±)p2,b{q±)) 



x 



P±(ki 



(45) 



^2,1 



^-^+ ( 5P-^+( 27r ) 4 5(p+)5(A;- -p-)5(g+)5(A;- - g~ 



(k ± - P± nk A 



p\q\ 



(46) 



S»-S»+SP+S°-(27i) 4 5(p + )5(k- -p-)5(q-)5(k + - g H 
(pL(fcl -P±)Pifi(Q±)) (p{ a (P±)P2,b(k ± - q ± )) 



X- 



ql(k±-p±) 



P 2 ±(k± - q±) 



(47) 



The other possibilities like Ai 12,2 an d A.2 2,1,1 do not contribute because they give terms 
proportional to 5(/c ± )5 2 (A;_|_). Because /c 2 = M 2 , these delta functions have no support. From 
these expressions, we see clearly that the cross-section is related to averages on sources like 
Eq. (liZTj) . The averages can then be evaluated using any model available. However, it is 
convenient for phenomenological applications to first relate averages to unintegrated gluon 
distribution functions. 

These averages can be related to the unintegrated distribution function of /c_i_-factorization 



like 



47 



48] 



(p\ a (p±)pi,b(q±)) 



4n 2 5, 



ab 



(N? - 1) 



P± + q± 



x / d 2 x ± e l(p± - q ^ x± yU 2 1 ' 



d 2 x± 



(48) 



4Ti 2 5 ab 
(N? - 1) 

x / d 2 y±e 



P± + q± 



J( P ±-g±Y(y±+b ± ) d( j )2 ( P± 2 9± I^-l) 



(49) 



where ^^^^^ are the unintegrated gluon distribution functions per unit area. By defini- 
tion, they are related to the unintegrated distribution functions by 

0i,2 ip±) = I d' 2 y ± - 



2 d(j) li2 (p±\y±) 



d 2 y± 



(50) 



where the integration is on the transverse extent of the nuclei. 

Combining all these results together, it is possible to compute the cross-section. The 
integration on the impact factor J d 2 bj_ in the definition of the cross-section gives a delta 
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function (2ir) 2 6 2 (p± — q±) in the second source average Eq. 
straightforward calculation, we find that the cross-section is: 

da pp ^ hX 



SJ). After a lengthy but 



{2^f2E k 



d 3 k 



16vr 4 



P^uapik) f d 2 q±d 2 p± 



1)k 2 



(2vr) 4 



(2iry8 2 (p ± + q ± -k A 



H^( P± ,q ± )Hf(p ± ,q ± ) 

p±q± 



(51) 



This has exactly the same structure as the cross-section obtained in fcj_-factorization shown 
in Eq. (ITTj) but now, it is derived as the low density limit of CGC The longitudinal momen- 
tum fraction x dependence of the unintegrated distribution functions is introduced through 
quantum evolution. The weight functions [pi.,2] obeys a non-linear evolution equation in 
x, the JIMWLK equation [4( 



43J. This makes the unintegrated gluon distribution 



41 



42 



function x-dependent 0(gi,/i 2 ) — > 4>(x, q±, fi 2 ) and Eqs. (I5TI) and (fT7|) can be compared. 



V. PHENOMENOLOGY FOR PROTON-PROTON COLLISIONS 

Eq. fl5Tl) (or (jTTJ) ) can be used to study the phenomenology of f 2 production and to study 
different parameterizations of the unintegrated gluon distribution function. By contracting 
the indices of P^ va p(k)H 1 ^ \p±, q±)H"^(p±, q±), the cross-section at mid-rapidity is given by 



d 2 k±dy 



d 2 q ± d 2 p±(p 1 (x', p 2 ± , /i 2 )(f)2(x', q 2 ± ,fi 2 ) 



2n(N 2 -1)k 2 
x5 2 (k ± -p ± - q ± ) 

fcj , [P±(g± • k±) + q±{p ± ■ k±)] 
M 2 3M 4 p 2 ±ql 



1 + 777 + o„^2 „ 2 f ( 52 ) 



where x' = y — ^ k± . As shown in this expression, the precise measurement of the f 2 differ- 
ential cross-section is a direct probe of unintegrated distribution functions. The unintegrated 
distribution functions obey evolution equations like the BFKL or the CCFM equation. De- 
pending on the physical model used or the approximations involved in the solution of these 
equations, the unintegrated distribution function can be parametrized in various ways. To 
get an estimate of f 2 production that can be compared with experimental data at RHIC, 
we use standard parameterizations [9J 1 . The measurement of f 2 production can be used to 
put constraints on models of the unintegrated gluon distribution function by comparing the 
predictions of many approaches. In this paper, we use the following parametrizations. 



1 We would like to thanks H. Jung for handing us his FORTRAN routine CAUNIGLU which evaluates nu- 
merically all of these parameterizations. It can be found at |http: / / www.desy.de~.jung/ cascade/updf.html| 
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A. DIG (Derivative of the Integrated Gluon distribution function) 



By ignoring that the unintegrated distribution function depends on a factorization scale 
/i, it is possible to invert Eq. (jSj) to get 



(/)dig(x,P±) 



2 dxg(x,n 



2\ 



d/j, 2 

where xg(x,[i 2 ) is the gluon parton distribution function. 



(53) 



9 2 

£i —Pi 



B. JB (J. Blumlein) 

The JB parameterization is based on a perturbative solution of the BFKL equation 0, l^] . 
In this parameterization, the unintegrated gluon distribution function is written as 

JB (x,pi,/i 2 ) = / dzg(z,p ± ,/j 2 )-G (-,/j, 2 ) (54) 

Jx Z \Z / 

where G(x, /i 2 ) is the collinear gluon distribution function and 



^(^pi,/! 2 ) = <( — -/ (55) 

-iofya.ln(i)ln(^)J, ifpi>/! 2 



2PT 



In this last equation, J and Jo are the Bessel's and modified Bessel's functions of the first 
kind respectively and 6t s = 3a s /ii where a s is the strong coupling constant (a s ~ 0.25). 

C. CCFM (Catani,Ciafaloni,Fiorani,Marchesini) 

The unintegrated gluon distribution function have been calculated numerically by solving 
the CCFM 



51 



52 



3, |53(] evolution equations using a Monte Carlo method 54J. The initial 
conditions of the evolution equation are determined from a fit of the proton structure function 
F2(x,Q 2 ). The result of this procedure for many data sets is implemented in the routine 
CAUNIGLU written by H. Jung. 

D. KMR (Kimber, Martin, Ryskin) 

The KMR parameterization starts at a certain scale given by p]_ ~ 1 GeV 2 . The non- 
perturbative part is given by the MRST collinear distribution function. The unintegrated 
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5 10 15 20 25 
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FIG. 2: Numerical results of the differential cross-section as a function of transverse momentum 
\k± \ for different parameterizations of the unintegrated distribution function at RHIC energy (y/s = 
200 GeV). The results are shown for midrapidity (y = 0). 

distribution function is given by 0, Q 

QKMU^M ,r , > , ( 56 ) 
II P± > P_L 

where the function Gkmr{x,P±, now depends on the scale /i and can be evaluated nu- 
merically. 

For numerical computations, we set the number of color to iV c = 3, the center of mass 
energy to y/s ~ 200 GeV (RHIC) and the mass of to 1.27 GeV. For the DIG parametriza- 
tion, we use the GRV NLN collinear gluon distribution function and for the JB parametriza- 
tion, we use the MRS collinear distribution. The numerical results of the differential cross- 
section at midrapidity for seven different parameterizations at RHIC energy are presented 
in Fig. [2j In Fig. [3j we present the results for the fc_i_-integrated (for < < y/s) 
differential cross-section for the CCFM parameterizations. All the numerical integrations 
are done with the CUBA package [621 using both CUHRE and VEGAS algorithms. 

There are some qualitative differences between the predictions of the different uninte- 
grated distribution functions. As shown in Fig. [21 the CCFM parameterizations are de- 
creasing much faster at large \k±\ compared to other parameterizations. At low \k±\, the 
shape of the curves is different, especially the parameterization CCFM B0. These features 
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FIG. 3: Numerical results of the differential cross-section as a function of rapidity y for the CCFM 
parameterizations of the unintegrated distribution function at RHIC energy (^/s = 200 GeV). The 
transverse momentum k± has been integrated. 

give us a way to discriminate between the unintegrated distribution functions. As seen in 
Fig. [3l the parameterizations also lead to results that have different magnitudes. However, 
there is still an uncertainty on the overall magnitude because of the value of the coupling 
constant that have to be fixed by experiments. We now discuss this crucial point. 

It is possible to compare our predictions to experimental data to fix the value of coupling 
constant k. The STAR collaboration have measured the production of f'2 mesons in the 
invariant mass spectrum of pions in proton-proton collisions at y/s ~ 200 GeV and midra- 



pidity 37|. The present data only show the invariant mass distribution in the transverse 



momentum bin 0.6 GeV < \k±\ < 0.8 GeV. From the graph presented in 37J, it is possible 



to extract a rough estimate of the average number of fi produced per collision. Taking into 
account that the branching ratio of — > titt is BRf^^ ps 0.85 20j, the average number 



of /2 produced per collision is n = f® 05 dy Jqq G ^ v d 2 k± d2 ^ dy ~ 0.0024. The total inelastic 
cross-section for pp-collisions at ^/s = 200 GeV is cr inc i rs 40 mb. Thus, we can convert the 
average multiplicity to a cross-section and we get that 

r0.5 r0.8 GeV 

a ps Oinei / dy I (j 2 fcj_ p^ 0.095 mb (Experimental) (57) 

i-0.5 io.6 GeV d k±ay 

This is an estimation of the experimental result for the f% production cross-section. 

We can then compute the integrated cross-section a using Eq. (TIT]) . The results are 
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Parameterization 


K 2 a 


K (GeV) 


CCFM JS2001 


0.0767 


0.558 


CCFM J2003 set 3 


0.142 


0.759 


CCFM AO 


0.134 


0.737 


CCFM BO 


0.178 


0.850 


DIG 


0.0424 


0.415 


JB 


0.0578 


0.484 


KMR 


0.0708 


0.536 



TABLE I: Numerical result for the total integrated cross-section using different parameteriza- 
tion of the unintegrated distribution functions. To compare with STAR data, we integrated over 
0.6 GeV < \k±\ < 0.8 GeV and —0.5 < y < 0.5. From this result, it is possible to fix the value of 
the coupling constant k. 



presented in Table HI This can then be used to determine the value of the coupling constant 
k and we find that depending on the parametrization used, k is between 0.415 GeV and 
0.850 GeV. This is of the same order of magnitude but much larger that k ~ 0.1 GeV 
obtained from pion decay by assuming a universal coupling of f 2 to the energy momentum 



tensor of all strongly interacting particles 3J, |36| or by tensor meson dominance 30J. The 



coupling of f 2 to gluons is much weaker than the coupling to pions. There are two reasons 
that explain this discrepancy. First, we are assuming that the f 2 meson can be produced 
as a stable on-shell particle. In f 2 decay, it is obviously not the case, the f 2 is unstable and 
decays in pions (or any decay product). Second, there are no justification a priori that the 
partonic sector should have the same coupling constant as that of the hadronic sector. In 
that sense, it is expected that the coupling of f 2 to gluons and hadrons will be different. 
Once the coupling constant is determined, all the curves presented in Fig. [2] and [3] become 
predictions and can be compared to experimental data. 



VI. CONCLUSION 

Using our effective theory based on a coupling between the f 2 meson and the energy- 
momentum tensor, we were able to compute the inclusive cross-section of f 2 tensor mesons 
in proton-proton collisions. This was done in the /c_|_-factorization formalism and in the 
CGC formalism. We showed that in the low density limit where the saturation effects are 
not taken into account, the result of the calculation of the cross-section in the CGC formalism 
reduces to the result of fcj_-factorization. In some sense, we use the CGC to "prove" Eq. 
for f 2 production. A rigorous treatment of factorization would involve a separation 
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between the long and short distance dynamics using resummation and a reorganization of 
the perturbative expansion at all orders. Of course, this is a hard task to perform. In the 
CGC, the proof reduces solely to a power counting in the color charge densities. This kind of 



result was expected as it was shown in other explicit calculations, namely for gluon 18l. Il9] 
and heavy quark production [171 ]. 

We also looked into the phenomenology of /2-meson production in proton-proton collisions 
at RHIC. We computed the differential cross-sections for many different parameterizations of 
the unintegrated distribution function. We used STAR data to fix the coupling constant and 
found that it is larger than in other approaches, even if it has the same order of magnitude. 
These numerical results, although they should be supplemented by more experimental data, 
can be used for the determination of the most accurate unintegrated distribution function. 
This is very interesting since /c^-factorization is one of the main computationnal tool for 
particle production (like heavy quarks and Higgs bosons [sj]) in high-energy 

collisions at RHIC,LHC and Tevatron. Having a consistent and precise unintegrated gluon 
distribution function is a crucial element in this formalism. Up to now, the main processes 
used to determine this distribution function were heavy quarks, dijets and gauge bosons pro- 
duction in pp and electron-proton collisions 61|. Many improvements have been performed 
to fit experimental data from these distributions, but there are still many uncertainties, 
especially at small k±. This uncertainty can also be seen in the variability of our numerical 
results for different parameterizations since they essentially lead to very different predictions. 
The /2 production gives another observable that can be used to put constraints on the value 
of unintegrated distribution functions. 

Our work is also a starting point for further studies such as the production of /2 in proton- 
nucleus collisions. According to CGC ideas, saturation effects should play an important role 
in that kind of process and thus, the /2 could be used to look at initial state effects. In that 
sense, our calculation can serve as a point of comparison for the presence and the magnitude 
of these initial state effects. This is under investigation and should be the topics of a future 
publication. 
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APPENDIX A: FEYNMAN RULES 



In this Appendix, the Feynman rules for gluon-/ 2 interaction are presented. By referring 
to the interaction Lagrangian Eqs. ([2]),(jl]and (j3j), it is easily seen that there are three possible 
interaction vertices: gg — > f 2 , ggg —>■ f 2 and gggg — > f 2 . However, the second and third 
ones are 0(g) and 0(g 2 ) respectively. They are not part of the leading order contribution, 
so they are not considered in this analysis. The vertex gg — ► f 2 can be evaluated and is 
given by 



K 



ab 



g fMT g up (q- P )+g pp g^(q-p) 



-g pu g ap {q ■ p) - g pa q p p v - g pp q v p a - g up q p p a - g va q p p 
+g pa q p p u + g pu q v p p + g pu q p p a 



(Al) 



where k ~ 0(100 MeV) to be fixed by experiments. The vertex has the following important 
property: 

P P V a 7 Pa (k,p, q) = q,V P r(k,p, q) = (A2) 

which is related to the conservation of the Abelian part of the energy-momentum tensor. 
The external f 2 have the polarization e^(fc). The f 2 Feynman propagator is 



where P^ pa (p) is the projection operator given by Eq. ( 113 



APPENDIX B: COLLINEAR FACTORIZED CROSS-SECTION 

The tensor meson production from pp-collision is investigated in this appendix, using 
the usual parton model (the collinear factorization formalism). This formalism can be used 
when Aq CD <C fi 2 ~ s where again s is the squared center of mass energy, Aqcd ~ 200 MeV 
is the usual QCD scale and fi is the typical parton interaction scale. At RHIC or LHC 
energy, this inequality is recovered at very large transverse momentum, when /i 2 « M\ ~ s, 
which may not be physically observable. However, it is interesting to make the calculation 
using this method as a consistency check for the fc^-factorized formalism. 

The cross-section for / 2 - m eson production in the collinear factorization formalism is given 

by 

(2n) 3 2E k = J dx 1 dx 2 G p \x 1 ,Q 2 )G p %x 2 ,Q 2 ) 

X{27T)32Ek ^k~ (B1) 
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where G P1 ' P2 are the gluon distribution function of the two protons, x\ 2 are the longitudinal 

p + 

momentum fraction (x = + art °" ) and Q is the factorization scale. The cross-section for 

Phadron 

on- shell gluons to on-shell tensor meson (gg — > f 2 ) is 

^ f2Ek ^k~ = ^\M 99 - h \ 2 {^\k-p-q) (B2) 

where Ai 99 ^? 2 is the matrix element for the process gg — > f'2 (see figured] and take on-shell 
gluons). The matrix element can be easily computed to lowest order from the Feynman 
rules described in Appendix |A] 

The matrix element squared is averaged over the degrees of freedom of the initial state 
and summed over the degrees of freedom of the final state. Once conservation of energy and 
momentum are used, it is given by 

- j^y E I E I^P = <B3) 

K c > a,b A,A',A" V c ' 

where A, A' and A" are the polarizations of the gluons and /2" m esons. This is obtained using 
the fact that for on-shell gluons, the sum on polarization is 

e,(k)e u (k) = - 9iXV + — + j-fLp - - 9yv (B4) 

A 

The replacement of the sum over polarizations by g^ (as depicted in the second part of 
the equation) is possible only because of the property (IA2[) of the vertex. The sum over 
polarizations of /2-meson is given by Eq. ( JT3l ). 

The initial gluons do not have transverse momentum in the collinear factorization for- 
malism (p± = q± = 0). Then, the momentum fractions are given by X\ = ^= and X2 = ^§ 
and we can write the four-momenta as 

P^^O^IP^l,-^ (B5) 



Using the expression for the matrix element and integrating the longitudinal components 
with the delta functions, the cross-section becomes 

(27r) 3 2ff fc ^ = _ (2rr) 2 6 2 (k ± )G(x + , Q 2 )G(x^ Q 2 ) (B6) 

The exact same equation was obtained by computing the collinear limit of fcj_-factorization, 
Eq. (jl~9l) . Similar results were obtained for rf- meson production 
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